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Abstract 

Small-amplitude weakly coupled oscillators of the Klein-Gordon lattices are approximated 
by equations of the discrete nonlinear Schrodinger type. We show how to justify this ap¬ 
proximation by two methods, which have been very popular in the recent literature. The 
hrst method relies on a priori energy estimates and multi-scale decompositions. The second 
method is based on a resonant normal form theorem. We show that although the two meth¬ 
ods are different in the implementation, they produce equivalent results as the end product. 
We also discuss applications of the discrete nonlinear Schrodinger equation in the context of 
existence and stability of breathers of the Klein-Gordon lattice. 


1 Introduction 

We consider the one-dimensional discrete Klein-Gordon (dKG) equation with the hard quartic 
potential in the form 

Xj -b Xj + x^ = e(xj+i — 2xj + Xj-i), j e Z, (1) 

where t E M is the evolution time, Xj{t) E M is the horizontal displacement of the j-th particle in 
the one-dimensional chain, and e > 0 is the coupling constant of the linear interaction between 
neighboring particles. The dKG equation ([1]) is associated with the conserved-in-time energy 

H = ]^'^x] + x) + e{xj+i - Xjf + ^ X 

which is also the Hamiltonian function of the dKG equation ([T]) written in the canonical variables 
{xj,Xj}j£i- The initial-value problem for the dKG equation ([T]) is globally well-posed in the 
sequence space ^^(Z), thanks to the coercivity of the energy in ([2|) in £^(Z). 

By using a scaling transformation 

Xj{i) = {I+ 2e)~^/‘^Xj{t), i = {I+ 2e)^/‘^t, e = (1-b 2e)“^e, (3) 
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and dropping the tilde notations, the dKG equation ([T]) can be rewritten without the diagonal 
terms in the discrete Laplacian operator, 

Xj + Xj + x^ = e{xj+i + Xj-i), j G Z. (4) 

Note that the values of e in Q are now restricted to the range (O, ^), because the map e —)• 
(1 + 2e)“^e is a diffeomorphism from (0,oo) to (O, ^). This restriction does not represent a 
limitation if we study solutions of the dKG equation for sufficiently small values of e. 

We consider the Gauchy problem for the dKG equation (|3]) and we aim at giving an approx¬ 
imation of its solutions by means of equations of the discrete nonlinear Schrodinger type, up to 
suitable time scales. This approach can be useful in general, but it may have additional interest 
when particular classes of solutions of the dKG equation Q are taken into account. In the case 
of systems of weakly coupled oscillators, relevant objects are given by time-periodic and spatially 
localized solutions called breathers. 

Existence and stability of breathers have been studied in the dKG equation in many recent 
works. In particular, exploring the limit of weak coupling between the nonlinear oscillators, 
existence |26) and stability of the fundamental (single-site) breathers were established (see 

also the recent works in mm)- More complicated multi-breathers were classified from the point 
of their spectral stability in the recent works [1] [25l [33]. Nonlinear stability and instability of 
multi-site breathers were recently studied in ca¬ 
lf the oscillators have small amplitudes in addition to being weakly coupled, the stability of 
multi-breathers in the dKG equation is related to the stability of multi-solitons in the discrete 
nonlinear Schrodinger (dNLS) equation: 

2iaj -|- 3\aj\‘^aj = -|- aj-i, j G Z, (5) 

where aj{et) G C is the envelope amplitude for the linear harmonic e** supported by the linear 
dKG equation (|1|) with e = 0. The relation between the dKG and dNLS equations (|1|) and ([5]) 
was observed in |29| based on numerical simulations and was elaborated in |33| with perturbation 
technique. 

The present contribution addresses the justification of the dNLS equation ([5|), and its gener¬ 
alizations, for the weakly coupled small-amplitude oscillators of the dKG equation (|4|). In fact, 
we are going to explore two alternative but complementary points of view on the justification 
process, which enables us to establish rigorous bounds on the error terms over the time scale, 
during which the dynamics of the dNLS equation ([5|) is observed. 

The first method in the justification of the dNLS equation ([5|) for small-amplitude weakly 
coupled oscillators of the dKG equation Q is based on a priori energy estimates and elementary 
continuation arguments. This method was used in the derivation of the dNLS equation [9| and the 
Korteweg-de Vries equation [a da dsiES] in a similar context of the Fermi-Pasta-Ulam lattice. 
The energy method is based on the decomposition of the solution into the leading-order multi-scale 
approximation and the error term. The error term is controlled by integrating the dKG equation 
with a small residual term over relevant time scale. The energy method is computationally efficient 
and simple enough for most practical applications. 

The second method is based on the resonant normal form theorem, which transforms the 
given Hamiltonian of the dKG equation to a simpler form by means of near-identity canonical 
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transformations BM- The normal form, once it is obtained in the sense of an abstract theorem, 
does not require any additional work for the derivation and justification of both the dNLS equation 
and its generalizations, which appear immediately in the corresponding relevant regimes. Starting 
from the works [mils], the normal form approach for the dKG equation was recently elaborated 
in |30| and applied in |31) for a stability result. 

We hope that the present discussion of the two equivalent methods can motivate readers for 
the choice of a suitable analytical technique in the justification analysis of similar problems of 
lattice dynamics. It is our understanding that the two methods are equivalent with respect to the 
results (error estimates, time scales) but have some differences in the way one proves such results. 

Besides justifying the dNLS equation ([5]) on the time scale we also extend the error 

bounds on the longer time intervals of 0(| log(e)|e“^). Similar improvements were reported in 
various other contexts of the justification analysis |9l[22]. Within our context, we will show in the 
end of this paper how to transfer the known results on existence and stability of multi-solitons 
in the dNLS equation ([5]) to the approximate results on existence and stability of multi-breathers 
in the dKG equation (j4]). In particular, we will address the relevant approximations for the dKG 
breathers on the extended time intervals obtained from the analytical results on the asymptotic 
stability of dNLS solitons Ha [Ml EH], quasi-periodic localized solutions [HIET], internal modes 
of dNLS solitons |32| . and nonlinear instability of multi-site solitons |23) . 

We finish introduction with a review of related results. Small-amplitude breathers of the dKG 
and dNLS equations were approximated with the continuous nonlinear Schrodinger equation in 
the works [ZlIHlEe]. An alternative derivation of the continuous nonlinear Schrodinger equation 
was discussed in the context of the Fermi-Pasta-Ulam lattice [HEQlEIlET]. In the opposite di¬ 
rection, derivation and justification of the dNLS equation from a continuous nonlinear Schrodinger 
equation with a periodic potential were developed in the works |34[ I35| . Finally, justification of 
the popular variational approximation for multi-solitons of the dNLS equation in the limit of weak 
coupling between the nonlinear oscillators is reported in m- 

The paper is organized as follows. Section 2 reports the justification results obtained from the 
energy method and multi-scale expansions. Section 3 reports the justification results obtained 
from the normal form theorem. Section 4 discusses applications of these results for existence and 
stability of breathers in the dKG equation. 

2 Justification of the dNLS equation with the energy method 

In what follows, we consider the limit of weak coupling between the nonlinear oscillators, where 
e is a small positive parameter. We also consider the small-amplitude oscillations starting with 
small-amplitude initial data. Hence, we use the scaling transformation xj = where p 

is another small positive parameter. Incorporating both small parameters, we rewrite the dKG 
equation (|3|) in the equivalent form 

'^3 + ^3 + ~ ^ ( 6 ) 

The standard approximation of multi-breathers in the dKG equation ([6|) with multi-solitons 
of the dNLS equation ([5|) corresponds to the balance p = £■ In Sections 2.1-2.3, we generalize the 
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standard dNLS approximation by assuming that <C p < e. In Section 2.4, we discuss further 
generalizations when p belongs to the asymptotic range <C p < e^. 


2.1 Preliminary estimates 

To recall the standard dNLS approximation, we define the slowly varying approximate solution 
of the dKG equation ([6]) in the form 

Xj{t) = aj{et)e^^ + aj{et)e~^^. (7) 

Substituting the leading-order solution ([7]) to the dKG equation ([6]) and removing the resonant 
terms at the leading order of 0{e), we obtain the dNLS equation in the form 

2iaj + 3u\aj\'^aj = Oj+i + Uj-i, j G Z, (8) 


where the dot denotes the derivative with respect to the slow time t = et and the parameter 
V = pjeis defined in the asymptotic range e <^v < 1. 

With the account of the dNLS equation ([8]), the leading-order solution ([7]) substituted into 
the dKG equation ([6]) produces the residual terms in the form 

Resj(t) := p (oje®* -|- aje“®*) . (9) 

The second residual term is resonant but occurs in the higher order which is not an obstacle 

in the justification analysis. The hrst residual term is non-resonant but it occurs at the leading 
order of 0{p) S> Therefore, the hrst term needs to be removed, which is achieved with the 

standard near-identity transformation. Namely, we extend the leading-order approximation ([7]) 
to the form 

Xj{t) = aj{et)e^^ + aj{et)e~'^^ H—p (a^(et)e^®* -|- o^(et)e“^®*) . (10) 

8 ^ j j ■> 

For simplicity, we do not mention that Xj depends on e and p. Substituting the approximation 
(11011 into the dKG equation ([6]), we obtain the new residual terms in the form 


Resj(t) := (dj-e®* -b a^e ®*) - ^ep ((a^+i -b -b (a^+i + a^_i)e 
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-b^p^ (a^e®* -b aje ®*)^ (a|e^®* -b o|e ^®*) -b ^ep [i 


2 ^ 


2 ■ 

lUjUje — lajttje 


-b^p^ (oje®* -b aje ®^) (a^e^®* -b a|e ^®‘)^ -b ^e^p + a?-e ^®‘ 
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( 11 ) 


Note that all time derivatives of aj in the residual term m can be eliminated from the dNLS 
equation ([8|) provided that {ojljgz is a twice differentiable sequence with respect to time. For 
all purposes we need, it is sufficient to consider the sequence space ^^(Z). Hence we denote the 
sequence {aj}j£z in by a- 

The next results give preliminary estimates on global solutions of the dNLS equation ([8]), the 
leading-order approximation (jlOjl . and the residual term (jlljl . 
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Lemma 1 For every ag € £^(Z) and every S M, there exists a unique global solution a{t) of the 
dNLS equation ^ in for every t E M such that a(0) = ao- Moreover, the solution a[t) is 

smooth in t and ||a(t )||^2 = ||ao||£ 2 . 

Proof. Local well-posedness and smoothness of the local solution a with respect to time variable 
t follow from the contraction principle applied to an integral version of the dNLS equation ([8|). 
The contraction principle can be applied because the discrete Laplacian operator is a bounded 
operator on whereas is a Banach algebra with respect to pointwise multiplication and 

the ^^(Z) norm is an upper bound for the ^°°(Z) norm of a sequence. Global continuation of the 
local solution a follows from the £^(Z) conservation of the dNLS equation dH). ■ 

Lemma 2 For every ao G there exists a positive constant C'x(||ao||^ 2 ) (that depends on 

||ao||£ 2 j such that for every p G (0,1] and every t G M, the leading-order approximation HI0(1 is 
estimated by 

||X(t )||,2 + ||X(t )||,2 <Cx(||ao||, 2 ). (12) 


Proof. The result follows from the Banach algebra property of l'^(Z) and the global existence 
result of Lemma [TJ ■ 

Lemma 3 Assume that p < e. For every ag G ^^(Z), there exists a positive e-independent constant 
C/jdlaoll^z) (that depends on ||ao||£ 2 j such that for every e G (0,1] and every t G M, the residual 
term in m is estimated by 

||Res(t )||^2 <CR{\\ao\\p)e'^. (13) 


Proof. The result follows from the Banach algebra property of ^^(Z), as well as from the global 
existence and smoothness of the solution a(t) of the dNLS equation ([8]) in Lemma [TJ ■ 

2.2 Justification of the dNLS equation on the dNLS time scale 

The main result of this section is the following justification theorem. 

Theorem 1 Assume that p is defined in the asymptotic range <C p < e. For every tq > 0, 
there is a small eg > 0 and positive constants Cq and C such that for every e G (0, eo), for whieh 
the initial data satisfies 

lle(O) - x(0)||p + \\m - X(0)|b2 < Cop-h^ (14) 

the solution of the dKG equation ^ satisfies for every t G [—Top~^,Top~^], 

lle(t) - X(t)|b2 + \m - X(t)|b2 < Gp-V. (15) 
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Remark 1 If p = e, the justification result of Theorem{l\ guarantees that the dynamics of small- 
amplitude oscillators follows closely the dynamics of the dNLS equation on the dNLS time 
scale [—To, To] for the variable r = et. 

Remark 2 If p = e®/®, the error term in I115\) satisfies the bound. The error term is 

controlled on the longer time scale [—roe”®/®] for the variable t = et of the dNLS equation 
with V = e®/®. 

To develop the justification analysis, we write f{t) = X(t) + y(t), where X(t) is the leading- 
order approximation (IlOp and y(t) is the error term. Substituting the decomposition into the 
lattice equation ([6]), we obtain the evolution problem for the error term: 

Vj + Vj + + ^p^jVj + pvj - 4yj+i + Vj-i) + Resj = o, j e z, (le) 

where the residual term Res(t) is given by (jlip if a(t) satisfies the dNLS equation ([8]). Associated 
with the evolution equation (I16p . we also define the energy of the error term 

[yl + yj + - ‘^^yjyj+^l ■ 


For every e G (O, |), the energy E{t) is coercive and controls the £"^{1,) norm of the solution in 

l|y(/)ll.^2 + ||y(t)||,^. <4F;(t), (18) 

for every t, for which the solution y{t) is defined. The rate of change for the energy (1171) is found 
from the evolution problem (I16p : 


dE 

dt 


^ [-yjReSj SpXjXjy'j - 3pXjy]yj - py^yj 
j&Z 


(19) 


Thanks to the coercivity (fTHl) . the Cauchy-Schwarz inequality, and the continuous embedding 
ofe{Z) to we obtain 


dE 

< 2E^/^ 

dt 



< [||Res(t) 11,2 +6p.F^/2||x(t) 11^, ||x(t) 11,2 + 12 pE||X(t )||,2 + 8pE®/2 


( 20 ) 


To simplify analysis, it is better to introduce the parametrization E = and rewrite (|2Up in the 
equivalent form 


dQ 

dt 


< ||Res(t )||,2 + 6 pQ||X(t)||, 2 ||X(t )||,2 + 12 pQ 2 ||x(f )||,2 + 8 pQ®. 


The energy estimate (1211) is the starting point for the proof of Theorem [TJ 


( 21 ) 


Proof of Theorem [H Let tq > 0 be fixed arbitrarily but independently of e and assume that 
the initial norm of the perturbation term satisfies the following bound 

Q(0) < Cop-h^ (22) 
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where Cq is a positive e-independent constant and e € (O, is sufficiently small. Note that the 
bound follows from the assumption ()14p and the energy (I17p subject to the choice of constant 
Co. 

To justify the dNLS equation ([8]) on the time scale for t, we define 

To := sup i fo € [0, ro/o“^] : sup Q{t) < CQp~^e^\ , (23) 

[ t6[-4o4o] J 

where Cq > Co is a positive e-independent constant to be determined below. By continuity of 
the solution in the norm, it is clear that Tq > 0. 

By using Lemmas [2] and [3l as well as the definition (I23D , we write the energy estimate (I2ip 
for every t S [—To,To] in the form 

< C/je^ -|- p (6C^ -|- VICxCqP + ^CqP ^e^) Q. (24) 

If e > 0 is sufficiently small and <C p, for every t G [—To,To], one can always find a positive 
e-independent ko such that 


dQ 

dt 


QCl + l2CxCQp-^e^ + 8Clp-‘^e^ < ko. (25) 

Integrating (1241) . we obtain 

g(t)e-^^o|t| _ g(o) < / Cne^e-P^^Ut'<=^. (26) 

Jo P^o 

By using ([2^ . we obtain for every t E [— ro,ro]: 

Qit) < p-h^ (Co + k^^Cn) (27) 

Hence, we can dehne Cq := (Co -|- and extend the time interval in (I23p by elementary 

continuation arguments to the full time span with Tq = tqp~^. This completes justification of the 
dNLS equation ([8|) in Theorem [T] □ 

2.3 Justification of the dNLS equation on the extended time scale 

Next, we justify the dNLS equation ([8]) on the extended time scale 

[-A\log{p)\p-\A\log{p)\p-^], (28) 

for the variable t, where the positive constant A is fixed independently of e. The main result of 

this section is the following justification theorem. 

Theorem 2 Assume that there is a ^ (0,1) such that p is defined in the asymptotic range 

2 

ei+“ <C p < e. 
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For every A £ (O, where /cq is defined in \33\) below, there is a small eo > 0 and positive 

constants Cq and C such that for every e £ (O,eo), for which the initial data satisfies 

lle(O) - X(0)|b2 + \\m - X(0)|b2 < Cop-h^ (29) 

the solution of the dKG equation ^ satisfies for every t in the time span i28\) . 

iie(t) - my + iie(t) - my < cp-m. m 


Remark 3 If p = e, the extended time scale h2^) corresponds to the interval [—^| log(e)|, log(e)|] 

for the variable t = et in the dNLS equation hence it extends to all times r as e ^ 0. 

Remark 4 If p = then the error term in LSfl\) satisfies the bound, which is 

small if a ^ (O, The error term is controlled on the longer time scale 

-Tol log(e)|e"^/^ro|log(e)|e"^/® 

for the variable t = et of the dNLS equation with v = e^!^. 

Proof. We use the same assumption (1221) on the initial norm of the perturbation term. To 
justify the dNLS equation ([8]) on the time scale (I28p for t, we define 

Tg := sup I to G [0,^1 log(p)|p"^] : sup Q{t) < CQp~^~°e‘^ I , (31) 

[ t&[-to,to] J 

where Cq is a positive e-independent constant to be determined below. 

By using Lemmas [2] and [3l as well as the definition (13ip , we write the energy estimate (I2ip 
for every t £ [—Tq,Tq] in the form 

< Crc^ + p (eCi + 12CxCQp-^-^e^ + 8C^p-2(i+“) e^) Q. (32) 

If e > 0 is sufficiently small and e^ <C p^~^°‘, then for every t £ [—TqjTq], one can always hnd a 
positive e-independent fco such that 

6Ci + 12CxCQp-i-“e2 + 8C^p-2(i+“)e^ < ko. (33) 

By integrating the energy estimate (f3^ in the same way as is done in (1261) . we obtain for every 
t £ [-To*,To*]: 

Q{t) < p-^e2(Co + A:o"^C/j)e'=“^l^°s(p)l 

< p-^-^e^ {Co + kfi^Cy , (34) 

where the last bound holds because koA £ (0,a). Hence, we can define Cq := Co + k^^Cn and 
extend the time interval in (I3ip by elementary continuation arguments to the full time span with 
Tq = A\ log(/o)|p“^. This completes justification of the dNLS equation (|8]) on the time scale ([28l) . 


dQ 

dt 
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2.4 Approximations with the generalized dNLS equation 


Extensions of the justification analysis are definitely possible by including more e-dependent terms 
into the dNLS equation ([8]) and the leading-order approximation (|10p . which makes the residual 
term (HID to be as small as ©(e”) for any n > 2. These extensions are not so important if 
p < e but they become crucial to capture the correct balance between linear and nonlinear 
effects on the dynamics of small-amplitude oscillators if p < e^. 

To illustrate these extensions, we show how to modify the justification analysis in the asymp¬ 
totic range <C p < e^. We use the same leading-order approximation (llOp in the form 

Xj{t) = aj{et)e''^ + aj{et)e~^^ H— p (a^(ef)e^** -|- , (35) 

8 ^ J j ■> 

but assume that a(r) with t = et satisfy the generalized dNLS equation 

“licLj -|- = Oj+i T i T ^ (flj-i-2 T 2(1^ -|- (ij— 2 ) i 3 S (36) 


Here we have introduced the parameter 5 = p/e^ in the asymptotic range e <C 5 < 1. Substituting 
(I35p and (I36p into the dKG equation ([6|), we obtain modifications of the residual terms (|lip in 
the form 

1 ■ 1 .. 

ReSj(f) := — (4hj -|- nj-1-2 + -|- Oj— 2 ) (hj -|- Oj +2 + 2aj -|- O‘j—2) c ** 

((“j-i-i + ^*0 

+ -b -b -b ^ep 


-,3„3it I 7^3„—3it\‘^ I 

8' 


-b—p^ (oj-e** -b Oj-e **) -b a|e -b -e^p -b o^e 


,3„3it I ■773„—3it 


(»?«“ + 


(37) 


By using the extended dNLS equation (|36p . we realize that the residual terms of the Op{e'^) order 
are canceled and the residual term in ()37p enjoys the improved estimate 


Res(t )||^2 < C'/^(||ao||£ 2 )e^ 


(38) 


compared with the previous estimate (|13p . As a result, the justification analysis developed in the 
proof of Theorems [T] and [2] holds verbatim and results in the following theorems. 


Theorem 3 Assume that p is defined in the asymptotic range p < . For every tq > 0, 

there is a small Cq > 0 and positive constants Cq and C such that for every e G (0, Cq), for which 
the initial data satisfies 

lie(O) - X(0)||p + \\m - X(0)|b2 < Gop-'e^, (39) 

the solution of the dKG equation ^ satisfies for every t G [— Top~^, Top~^], 

\\m - x(f)ib2 + \m - x(f)ib2 < cp-h^ (40) 
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Theorem 4 Assume that there is a & (^> 5 ) that p is defined in the asymptotic range 

p < e^. 

There is ^0 > 0 such that for every A S (0, Aq), there is a small eo > 0 and positive constants Cq 
and C such that for every t G (0, € 0 ), for which the initial data satisfies 

11^(0) - X{0)\\p + ||<e(0) - x(0)||p < Cop-h^ (41) 

the solution of the dKG equation ^ satisfies for every t in the time span h2^) . 

Il^(i) - X(i)||p + \m - X(i)|b2 < Cp-^-^e\ (42) 


We note that X in Theorems [3] and [4] is defined by the leading-order approximation (f35l) . 
whereas a satisfies the generalized dNLS equation (I36p . The time scales in Theorems [3] and 0] are 
appropriate for the generalized dNLS equation (|36p because < 1 and €p~^ > 

3 Justification of the dNLS equation with the normal form method 

We now consider the dKG equation (01 on a finite chain of 2N + 1 oscillators under periodic 
boundary conditions. The finite dKG chain is associated with the Hamiltonian H = Hq + Hi, 
where 

1 iV ^ N 

Ho-.= - [y| + x] - 2exj+iXj\ , 4 X] 

3=-N j=-N 

subject to the periodic boundary conditions X-n = xtv+i and = yN+i- Since N can be 
considered arbitrary large in the perturbation approach which follows, the finite dKG chain ap¬ 
proximates the infinite problem in the asymptotic sense, as N —)• 00 . Although the present theory 
can be adapted to the infinite lattice, we prefer to rely on some already proved results for the 
finite dKG chain for the sake of brevity. 

According to the previous result in |30| . for any small coupling e, there exists a canonical 
transformation Tx which puts the Hamiltonian H = Hq + Hi, with Hq and Hi in (j43p . into an 
extensive resonant normal form of order r 

=Hn + Z + , {Hn, Z} = 0 , (44) 

where Hq is the Hamiltonian for the system of 2N + 1 identical oscillators of frequency H (which 
is the average of the linear frequencies |18)). Z is a non-homogeneous polynomial of order 2r -|- 2, 

p(r+l) 

is a remainder of order 2r -|- 4 and higher, and r grows as an inverse power of e. Such a 
normal form was shown to be well defined in a small ball H^i/2(0) C P of the phase space V, 
endowed with the Euclidean norm (which becomes the ^^(Z) norm in the limit N —)• 00 ), provided 
<C 1. The linear part of the Hamiltonian Hq = Hp is equivalent to the selected squared 
norm (uniformly with N), thus the almost invariance of Hq over times |t| ~ (r^/o)“'’“^ is easily 
derived since Hq = {Hq,P^^~^^'^}. 
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Looking at the structure of the normal form H^ + Z produces a generalized dNLS equation, 
where all oscillators are coupled to all neighbors and the coupling coefficients both for linear and 
nonlinear terms decay exponentially with the distance between sites. To be more specific, Z can 
be split as the sum of homogeneous polynomials Zq, Zi, where Zq is quadratic, Zi is quartic, 
and Zr is of the order 2r + 2 > 4. Each of these homogenous polynomials can be developed in 
powers of the coupling coefficient e, where the term of order e™ is responsible for the coupling 
between lattice sites separated by the distance m. The key ingredient to obtain the normal form 
is the preservation of the translation invariance (called cyclic symmetry in |181 130)). which also 
allows us to produce estimates that are uniform with N. 

If we limit to r = 1, the transformed Hamiltonian (|44|) reads 

H-(i) = /C + P(2) ^ ]C:= Hn + Zo + Zi , 

where the quadratic and quartic polynomials Zq and Zi include all-to-all interactions, exponen¬ 
tially decaying with e. Hence, K, represents Hamiltonian of the generalized dNLS equation. If we 
truncate both Zq and Zi at the leading order in e, we recover Hamiltonian of the usual dNLS 
equation. 

In Section 3.1, we introduce some definitions. The linear transformation is analyzed in Section 
3.2. The nonlinear normal form transformation is performed in Section 3.3. Approximations with 
the usual dNLS equation are obtained in Section [3.41 Approximations with the generalized dNLS 
equation are discussed in Section 3.5. 

3.1 Some definitions 

We start recalling some definitions which allow us to characterize the structure of the normal form 

(see also [13 dSl [301 EH). 

Cyclic symmetry: We formalize the translational invariance of the model (1431) by using the 
idea of cyclic symmetry. The cyclic permutation operator r is defined as 

t{x-n, ■ ■ ■, xn) = {x-N+1, • • •, xn,x-n)- (45) 

This operator can be applied separately to the variables x and y. We extend the action of this 
operator on the space of functions as {Tf){x,y) = f{TX,Ty). 

Definition 1 We say that a function F is cyclically symmetric if tF = F. 

We introduce now an operator, indicated by an upper index 0, acting on functions: given a 
function /, a new function F = /® is constructed as 

N 

F = f^:= • (46) 

l=-N 

We say that f®{x,y) is generated by the seed f(x,y). Our convention is to denote the cyclically 
symmetric functions by capital letters and their seeds by the corresponding lower case letters. 
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Polynomial norms: Since we are interested in homogeneous polynomials (due to the use of 
Taylor expansion), we introduce the following norms. 


Definition 2 Let f{x,y) = fj,kX^y^ be a homogeneous polynomial of degree s in x, y 

and F = /®. Given a positive radius R, we define the polynomial norm of f by 

II/IIr := R‘ E 1.^^! ■ («) 

Ii| + |fc|=s 

Correspondingly, the extensive norm of F is given by 

||i"||?^ = ll/lk, (48) 


Vector fields: Let F be an extensive Hamiltonian with seed /; we will make use of the notation 
Xp to indicate the associated Hamiltonian vector field JXF, with J given by the standard Poisson 
structure. The Hamiltonian vector field inherits, in a particular form, the cyclic symmetry: indeed 
it holds true (see IMIEI]) 

d^^F = T^d^,F, dy.F = T^dy,F, j = -N,...,N. (49) 

As a result, a possible (but not unique) choice for its seed turns out to be the couple [dy^F, —dx^F). 
This fact allows us to define in a reasonable and consistent way the following norm 


Xf 


R 






\R 


(50) 


Interaction range and centered alignment: Let us now consider monomials x^y^ in multi¬ 
index notations for {j,k). 


Definition 3 Given the exponents {j,k), we define the support S{x^y^) of the monomials x^y^ 
and the interaction distance £{x^y^) as follows: 

S{x^y’") = {I : jifiO or ki 0} , i{x^y^) = diam(5(x^/)) . (51) 


We want to stress that, differently from what has been developed in [13 US], it is possible to 
implement tha same construction by asking the seeds of all the functions to be centered aligned, 
according to the following definition m- 


Definition 4 Let F = /® be a cyclically symmetric functions, with f depending on 2N + 1 
variables, f = f{x-N, ..., xq, ... ,xn). The seed f is said centered aligned if it admits the decom¬ 
position 

N 

S(/(-))C[-m,...,m] . (52) 

7n=0 
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Exponential decay: In order to formalize and control the interaction range, we introduce 

Definition 5 The seed f of a function F is said to he of class T>{Cf, fi) if there exist two positive 
constants Cf and p, < 1 such that for any centered aligned component it holds 


f 


(m) 


< CftT 


m = 0,... ,N 


3.2 Linear transformation 


Let us focus on the harmonic part Hq of the Hamiltonian H. From ()^3l) . Hq can be written as 
the quadratic form 

1 1 

Ho{x,y) =-y-y +-Ax ■ X (53) 

where H is a circulant and symmetric matrix given by 

A :=I — e(r + r''') . (54) 

Here r = {jij) is the matrix representing the cyclic permutation ()45]) . i.e. with Tij = (mod 2 Ar+i) 

using the Kronecker’s delta notation. 

Proposition 1 For every e £ (0, ^) the canonical linear transformation q = A^^^x, p = A~^^^y 
transforms the quadratic Hamiltonian Hq to the quadratic normal form 


H'(o) =Hn + ZQ 


{Hq, Zq} — 0 , 


(55) 


where Hq = /i® and Zq = C,® are cyclically symmetric polynomials, with centered aligned seeds 
Kq and ^Q of the form 

(56) 

(57) 


hn = ^{Qq +Po) 


and 


N 


Co — ^ ^ Co \ Co ^ — bm[QoiQm A q—m) + PoiPm + P-m)]- 

7n=l 

Here and bm are defined by 

N+l 


n := 


1 


2N + 


j=-N 


; bfji .— 




(58) 


whereas coj are the frequencies of the normal modes of Hq. Moreover, there exists a suitable 
positive constant such that each component satisfies the exponential decay 


"M 


< Cco(2e)- , 


hence Co £ ■ 
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Proof. We give here only few ideas to grasp the exponential decay of the all-to-all interactions 
due to the linear transformation. After applying q = p = we have 

Ho = ^p'^A^/^p+^q'^A^/^q. (59) 

By defining T := r + r^, one can rewrite as 

A^/2 ^ (I _ g2^)i/2 ^ ^ 

1=0 ^ ' 

In order to obtain the decomposition (I55p , we separate the diagonal part from the off-diagonal part 
AAI‘^ = n H-B and insert this decomposition into (|59p . The exponential decay (2e)™' comes from 
the observation that = 0 for all 0 < Z < m and from the estimate \{T'^)i < 2™. 

One can restrict to consider only the first raw due to the circulant nature of all the matrices 
involved (for all details see Appendix 6.1.1 in |18)). ■ 

Proposition 2 Under the linear transformation in Proposition Q], the quartic part Hi given in 
is cyclically symmetric with a centered aligned seed Hi = hf given by 

N 

hi = Y, ■ ( 60 ) 

m=0 

Moreover, there exists a suitable positive constant Ch^ such that each component satisfies the 
exponential decay 

hence hi EP(C'/n,2e). 

We can translate Propositions [T] and [2] by saying that in a suitable set of coordinates, the 
coupling part of the quadratic Hamiltonian Ho shows all-to-all linear interactions, with an expo¬ 
nentially decaying strength with respect to the distance between the sites. Such a linear trans¬ 
formation introduces similar all-to-all interactions also in the quartic Hamiltonian Hi. Moreover, 
in the new coordinates qj, the seed hi of the quartic term has the same exponential decay as the 
seed ^0 of th® quadratic term does. 

3.3 First-order nonlinear normal form transformation 

With the Hamiltonian H transformed by means of Propositions [1] and [2] into the form 

H = Hu + Zo + Hi , (61) 

we are now ready to state the (first-order) normal form theorem. This first-order theorem rep¬ 
resents the easiest formulation of the more generic Theorem 1 of |30| . The idea is to perform, 
by using the Lie transform algorithm explained in m, one normalizing step, provided e is small 
enough. Moreover, the normalizing canonical transformation is well defined in a (small) neigh¬ 
borhood Bpi /2 of the origin, where p is sufficiently small. 


<ChA2er , 
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Theorem 5 Consider the Hamiltonian H = /i® + + /i® with seeds hfi, Co; (1561) . (1571) . 

and ini). There exist positive 7 , e* < | and C* such that for every e G (0, e*), there exists a 
generating function Ti = x® of a Lie transform such that = H, where is a cyclically 

symmetric function of the form 

77 ( 1 ) =Hn + Zo + Zi + P^^'> , (62) 

with 0 = {Hq, Zq} = {Hfi,Zi}, whereas Z\ = is a polynomial of degree four with a seed Ci 
is of class V (C'/i^, 2 e), and is a remainder that includes terms of degree egual or bigger than 
six. Moreover, if the smallness condition on the energy 


1 

— 96(1 + e)a 

is satisfied, then the following statements hold true: 


(63) 


1. defines an analytic canonical transformation on the domain ^2 1/2 such that 




C Tx^Bi T_/2 C B i/2 


Byi /2 c T^^Byi /2 C Bpi /2 


Moreover, the deformation of the domain B 2 1/2 is controlled by 

3 P 


z G B 2 1/2 

9.P 


\\Tx, {z) - z|| < 4 ^C'*p 3/2 , Tf,l {z)-z < 4 ^ap 3/2 . (64) 




the remainder is an analytic function on B 2 1 / 2 , and it is represented by a series of cyclically 

3 P 

symmetric homogeneous polynomials of degree 2 s + 2 


s=2 



GV{2Cr^Ch„VTe) . 


(65) 


The interval (0, e*) with e* < i comes from the inequality 


/ 311 yi-26)[l-(26)l 
[64C/;J VTe 


^ > 1 


^Co 

(see for reference formula (33) in |31|i. and the constants C* and 7 can be written as 


a = 


and 


7 = 211 1 - 


37(1 - 2e) 1 - (2e)4 

1 


2/(e) 


H < 7 < 211. 


( 66 ) 


(67) 


Since e is sufficiently smaller than i, the constants C* is essentially independent on e, i.e. 

'Ch, 


c, = o 


n 


which implies that the same holds true for the threshold p* so that 

211 


3Cftj (1 + e) 


( 68 ) 
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3.4 Approximation with the dNLS equation 

We apply here the normal form transformation of Theorem [5] in order to approximate the Cauchy 
problem i: = {H,z} of the finite dKG equation (j4]) with a small initial datum zq. Let us denote 
with K, := -ftfj + Zq + Z\ the normal form part of the Hamiltonian = /C + in formula 
(I62p . Since Zq and Z\ have centered aligned seeds with the exponential decay, see decompositions 
([57)1 and (l60]) . we have 


N 

Z„=Y. 4”’ . 4"*’ (cP) (69) 

?n=l 

and 

N 

Zr = J2 := ■ (70) 

m=0 

Note that the expansion for Zq starts at m = 1, while Zi starts with m = 0. By truncating the 
e expansion of each normal form term Zj at their leading orders, we define the effective normal 
form Hamiltonian /Cgff as 

/Ceff := Hn + 4'^ + , -^res := /C - /Cgfr . (71) 

As already stressed in m, the truncated normal form /Cgfr represents Hamiltonian of the dNLS 
equation. In complex coordinates ifj = [qj + ipj) 1^/2, Hamiltonian /Cgfr reads as 

/Ceff = (H + 26i) ^ - 6i ^ ^ X] 

3 3 3 

where 6 i = 0{e) <0 is the same as in the expression (I58p of Proposition [TJ The corresponding 
dNLS equation is 

= ^'^3 + + V'i-i) + , (73) 

oipj 4 

and it has the same structure as the dNLS equation ([ 8 ]). 

We denote with z{t) the evolution of the dKG transformed Hamiltonian /C + with Za{t) 
the evolution of the dNLS model /Ceff and consequently with 5{t) the error 

5{t) ;= z{t) - Za{t) . (74) 


The two time scales over which we control the error of the approximation are given by 


Tq: = 


1 

? 

P 



(75) 


where a € (0,1) is an arbitrary parameter, and kq = 0{Chf) is given in (|86p . Similar definitions 
are used in (|23p and fl3T]) , in the proof of Theorems [T] and [2j 
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Theorem 6 Let us take p fulfilling ([63|) and e € (0,e*) as in Theorem\^ Let us first consider 
the two independent parameters p and e in the regime <C /O < e. Then, there exists a positive 
constant C independent of p and e such that for any initial datum zq G i? 2 ^i /2 with ||(5o|| ^ > 

the following holds true 

m)\\<Cp-^ih\ |i|<ro. (76) 

2 

Let us now consider the two independent parameters p and e in the regime ei+“ p < e, where 
a G (0,1) is arbitrary. Then, there exists a positive constant C independent of p and e such that 
for any initial datum zq G - 62 ^ 1/2 with ||5o|| ^ following holds true 

||5(t)|| |t|<ro*. (77) 

Remark 5 The upper bound for the error 6 given in (I76D and (HZD refers to the time evolution of 
the normal form (I72p in the transformed variables ^|J, which are near-identity deformations of the 
original variables {x,y). Since the transformation Tx is Lipschitz, with a Lipschitz constant L of 
order L = 0(1), the same bound of the error holds also in the original coordinates. Thus, from the 
analytic point of view, the nonlinear deformation of the variables does not affect the dependence 
of the estimates on p and e: only the constant C is changed by the Lipschitz factor L. 

Remark 6 The above estimates are equivalent, both in terms of error smallness and time scale, 
to the ones obtained in Theorems\^ and\^ once the original variables Xj = are recovered. 

Remark 7 The requirement p on the time scale Tq is needed in order to provide a meaningful 
approximation, which means that the error is much smaller than the leading approximation Za(t) 

p(i)ll < ~ l|2a(i)ll • 

2 

The same reason lies behind the requirement ei+“ p on the extended time scale Tq . 

Proof of Theorem Following a standard approach (see a similar problem in [6]), we first 
decompose the Hamiltonian H = Hl + Rat in its quadratic and quartic parts 

so that /Ceff = Hl + Rat — — KL^es- Correspondingly, the vector field is decomposed as 

Xj{ = -\- Xjsf. Denote the linear operator for Xjj^ by C. The equation of motions for z{t) 

and Za(t) reads 

with Res(t) := Xp( 2 )(za(t)) + XK:,,fiza(t)) . (78) 

The error 6{t) defined by (1741) satisfies the equation 

6 = C6+ [X^iza + 5)- XN(za)] + Res(f) , (79) 


z = Cz + Xn{z) , 

Za = CZa + X]S[(Za) - Res(f) 
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whose solution, with the initial value Sq, is given by Duhamel formula 


6{t) = e^^So + [ e + 5) - XAr(za) + Res(s)]ds . (80) 

Jo 

Now, since = 0, one has that C is an isometry. This allows to estimate 

I|(^(i)ll < ll'^oll + f [||X 7 v(za(s) + (i(s))-XAr(za(s))|| + ||Res(s)||]ds . (81) 

Jo 


The second term in the r.h.s. can be estimated with the definition of the residual and using the 
information that Za{t) preserves the norm, as a consequence of the conservation of Hq 

II ^ /„^MI / r- ^ ^ ^ [Cc_oP^/‘^e^ + Chy/'^e] 

||-X^p(2) (2^a('5)) II — ^ ll'^^res ('^Ol(^)) II — C _ 26)^ ^ 

where the two contributions in the second inequality come from the truncation of Zq and Z\ 
respectively. Thus, we obtain 

a /2 


\\Re^[s)\\<Cj^^^^[Cy + Ch,pe + Cuyy] . 

On the other hand, if 

ll<^ll < Ikall ~ P^^‘^ , 

then the increment of the nonlinear field can be well approximated by 

\\XN{Za{s) + 5(5)) - X^iZaism < ||X)^(Ca)|| ||5|| , 

where 

Ca := + A5 , A G (0,1) . 

If the smallness condition fl83)l for 6 holds, then ||C(j|| ~ which implies 

||X^(^a(5) + 5(5)) - X^(Za(5))|| < ||X^||^ ||5|| . 
By using the decomposition X'j^ = + Xy~^ it is possible to obtain 

Ch, 


(82) 

(83) 


|^Af||pl /2 ^ ^ 1 " 


_ -p . (84) 

(1-^)2^ ^ ^ 

By inserting (j82jl and fl84)l into fl 8 T]l . one gets a typical Gronwall-like integral inequality (see, e.g.. 
Lemma A.2 in M), which provides the time-dependent upper bound 

||5(t)|| < ||5o|| + CpV2 [f! + e + c,p] (e"“^' - l) 

IP J 

< e^°P^p-^/‘^e^ + Cp -^/2 (gKopt _ + pg + c^p^] , (85) 

where kq provides an upper bound for in (|84p 

(1 _ ^^)2 = 

and C depends only on e*, C^p, Ch^. Then, the bound (I76h follows from the assumption p < e. 

The bound ()77p is obtained similarly, just replacing the time span Tq in the above ()85p . which 
easily provides the factor p~°‘ in front of the estimate. □ 
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3.5 Approximations with the generalized dNLS equation 

The standard dNLS approximation is no more valid when ~ p. Indeed, in such a case, the 
contribution e^p~^ coming from the truncation of the linear field in (|85p is of order one, 

hence the error 5{t) can be comparable with the approximation Za{t) 

||<5(t)|| < ~ lka(t)|| • 

( 2 ) 

In such a regime, it is then necessary to include in the Hamiltonian /Cgfr at least the term Zq , 
responsible for the next-nearest neighbourhood linear interaction: 

/Ceff := Hn + . (87) 

Following the same steps as in the proof of Theorem [ 6 l it is possible to prove the following result, 
which is fully equivalent to Theorems [3] and 01 

Theorem 7 Let us take p fulfilling ([ 6 ^ and e G (0,e*) as in Theorem\^ Let us first consider 
the two independent parameters p and e in the regime <C p < . Then, there exists a positive 

constant C independent on p and e such that for any initial datum zq G B 2 ^ with ||(5o|| < , 

it holds true 

||5(t)|| <Cp-vv, |t|<ro. ( 88 ) 

3 

Let us now consider the two independent parameters p and e in the regime ei+“ p < e, where 
a G (0,1) is arbitrary. Then, there exists a positive constant C independent of p and e such that 
for any initial datum zq G i ?2 1/2 with ||5o|| < the following holds true 

||<5(t)|| , |t|<ro*. (89) 

The result of Theorem [7] yields Hamiltonian for the generalized dNLS equation: 

/Ceff = (H + 2bi + 262 ) l^il^ - ^2 ^ l^i +2 + 

j j j j 

where 62 = C^(e^) < 0 is the same as in the expression (I58p of Proposition [T] The corresponding 
generalized dNLS equation is 

3 

iipj = Llfij + bi{fij+i -t- ifj-i) + h 2 {'ifj +2 + 'fij- 2 ) + > (91) 

which has the same structure as the generalized dNLS equation (|36l) . Indeed, remembering that 
H in (1911) also has an expansion in e, and that the time variable is rescaled with e in (1361) . we can 
rewrite the right-hand-side of the generalized dNLS equation (I36p as follows: 

-Oj + (oj + i + «i-i) + ^(oi +2 + aj- 2 ) ■ 

This shows an e correction to the nearest neighbour coefficient, which in the normal form approach 
is embedded in the e-dependence of H, 61 , 62 and of the transformed coordinates. 
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More generally, within the normal form approach, different regimes of parameters can be 
treated with no efforts: once the requested scaling between e and p is chosen, one easily derives the 
minimal, and also the optimal, number of terms in the expansions of Zq and Zi to be included. The 
estimates follows as easy as before. Here we give the estimates for a general choice of truncation: 




l-i 




n—l 


7 {3) 




j=o 


(92) 


where N > I >2 and N >n>l. The error term 5 is now estimated similarly to (f85|) as follows: 




+ e”" + C**/? 


(gKopt _ 


l>2 , 


n > 1 . 


(93) 


Hence one can deal with all the regimes and with the desired error precision in a compact and 
flexible way. The extension to higher order terms in the nonlinearity would require further steps 
of the normal form transformations, thus modifying thresholds e* and p*, following the general 
version of Theorem [5] given in |30| . 


4 Applications of the dNLS equation 

We conclude the paper with a brief account of possible applications of the dNLS equations dH]) and 
(17311 . and their generalizations (|36p and (I9ip . in the context of small-amplitude weakly coupled 
oscillators of the dKG equation (jl]). 

1. Existence of breathers. Breathers of the dKG equation (time-periodic solutions localized 
on the lattice) can be constructed approximately by looking at the discrete solitons of the dNLS 
equation ([5]) in the form aj{T) = where fl € M\[—1,1] is defined outside the spectral band 

of the linearized dNLS equation and A G .^^(M) is time-independent. 

The limit e —>• 0 is referred to as the anti-continuum limit of the dKG equation (j3|), when 
the breathers at a fixed energy are continued uniquely from the limiting configurations supported 
on few lattice sites |26l I33| . Gompared to the anti-continuum limit, the dNLS approximation is 
very different, because the discrete solitons of the dNLS equation ([5]) are not nearly compactly 
supported due to the fact that the dNLS equation (jS]) has no small parameter. In agreement 
with this picture, the continuation arguments in |26[ [33] are not valid in the small-amplitude 
approximation, when the breather period T is defined near the linear limit 27r, because the inverse 
linearized operators become unbounded in the linear oscillator limit as T —)• 27r. 

By Theorem [H discrete solitons of the dNLS equation m are continued as approximate 
breather solutions of the dKG equation (j4]), which are only periodic solutions up to the time scale 
0(e“^). However, this approximation can be extended to all times, by considering time-periodic 
solutions of the dKG equation (j4]) with small e, using Fourier series in time, and eliminating all but 
the first Fourier harmonic by a Lyapunov-Schmidt reduction procedure. Then, the components of 
the first Fourier component satisfies a stationary dNLS-type equation, where the dNLS equation 
(|5|) is the leading equation. In this way, similarly to the work |34| . one can justify continuation of 
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discrete solitons of the dNLS equation ([5]) as approximate solutions of the true breathers in the 
dKG equation (j4]). 

Within the same scheme of Lyapunov-Schmidt decomposition, another equivalent route to 
prove the existence of breather for the dKG equation (j4]) is obtained by means of Theorem [5l 
Indeed, the discrete solitons of the dNLS equation ([5]) can be characterized as constrained critical 
points of energy, which are continued, under non-degeneracy conditions, to critical points of the 
true energy of the dKG equation Q. 

2. Spectral stability of breathers. Spectral stability of breathers in the dKG equation (jH 
can be related to the spectral stability of solitons in the dNLS equation ([5]). By Theorem[Tl we are 
not able to relate stable or unstable eigenvalues of the dNLS solitons with the Floquet multipliers 
of the dKG breathers, because the error term also grows exponentially at the time scale 0(e“^) 
(the same problem is discussed in |14) in the context of stability of the travelling waves in FPU 
lattices). However, by Theorem [2] obtained on the extended time scale 0(e“^| log(e)|), we can 
conclude that all unstable eigenvalues of the dNLS solitons persist as unstable Floquet multipliers 
of the dKG breathers within the 0{e) distance from the unit circle. 

If unstable eigenvalues of the dNLS solitons do not exist, we only obtain approximate spectral 
stability of the dKG breathers, because unstable Floquet multipliers of the dKG breathers may 
still exist on the distance smaller the 0{e) to the unit circle. On the other hand, if the spectrally 
stable dNLS solitons are known to have internal modes |32], then the Floquet multipliers of the 
dKG breathers persist on the unit circle by known symmetries of the Floquet multipliers |33| . 

3. Long time stability of breathers. By means of the normal form approach, it is possible 
to prove the long time stability result for single-site (fundamental) breather solutions of the dKG 
equation (j3|). Indeed, the variational characterization of the existence problem for such breathers 
in the normal form essentially implies an orbital stability in the normal form, which is translated 
into a long time stability in the original dKG equation m- 

In the case of multi-site dNLS solitons, nonlinear instability is induced by isolated internal 
modes of negative Krein signature, which are coupled with the continuous spectrum by nonlinear¬ 
ity [23]. By using the extended time scale 0(e~^\ log(e)|) of Theorem [21 we can predict persistence 
of this instability for small-amplitude dKG breathers. This was recently confirmed for multi-site 
dKG breathers in m- 

Also quasi-periodic localized solutions were constructed for the dNLS equation in the situa¬ 
tion, when the internal mode of the dNLS soliton occurs on the other side of the spectral band 
of the continuous spectrum dlllST]. These solutions correspond approximately to quasi-periodic 
dKG breathers. It is still an open question to consider true quasi-periodic breather solutions of 
the dKG equation (|4|). 

4. Asymptotic stability of breathers. Asymptotic stability of dNLS solitons supported 

at a single-site linear potential was considered in within the dNLS equation with the 

seventh-order power nonlinearity. Again, these results apply to the dKG equation only approx¬ 
imately on the extended time scale 0(e“^| log(e)|) by Theorem [2j However, by using canonical 
transformations and dispersive decay estimates, asymptotic stability of single-site breathers was 
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proved recently in [1], also in the case of the seventh-order power nonlinearity. 

The multi-site dKG breathers with isolated internal modes of positive Krein signature, which 
are coupled to the spectral band of the continuous spectrum due to nonlinearity, are also ex¬ 
pected to remain nonlinearly (and, perhaps, asymptotically) stable |13| . It remains however an 
open problem to study this problem directly for the dKG equation ([1|), without rescuing to the 
approximation result of Theorem [2j 
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